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Spreading phenomena are ubiquitous in nature and society. For example, disease, rumor, and 
information spread over underlying social and information networks. It is well known that there is 
no threshold for epidemic models on scale-free networks; this suggests that disease can spread on 
such networks, regardless of how low the contact rate may be. In this paper, I consider six models 
with different contact and propagation mechanisms. Each model is analyzed by degree-based mean- 
field theory. I show that the presence or absence of an outbreak threshold depends on the contact 
and propagation mechanism. 


I. INTRODUCTION 

The development of communication technology and 
transportation has increased connectivity among peo¬ 
ple. Outbreaks of several new infectious diseases, such 
as AIDS, SARS, swine flu, and Ebola, have threatened 
human society. These diseases spread over networks 
of contacts between individuals. Computer viruses and 
worms spreading through the Internet have caused se¬ 
vere economic damages all over the world. Moreover, 
rumors, opinions and innovations spread through human 
networks. Understanding the intrinsic mechanism be¬ 
hind spreading phenomena in networks is an important 
and urgent taskf^-l^. 

Epidemic dynamics has often been described by ordi¬ 
nary differential equations, that assume that the proba¬ 
bility for an infected individual to encounter a susceptible 
host is uniform!^ Q. However, social networks are not 
uniformly mixed but are highly heterogeneous. Many so¬ 
cial networks, such as telephone calls Q, e-mailssex¬ 
ual relationships 0 , file actor coraboration[T^, citation 
networks (ll|. and the Internet [l^, have scale-free prop¬ 
erties. A network is called scale-free if the distribution 
of d egree ( i.e., the number of links that connect to a 
node |13l4l6l| ) obeys a power law: 

P{k) - 

where k represents degree. For most real world networks, 
the exponent 7 is between 2 and 3. High degree nodes are 
called hubs. Spreading processes in such networks have 
been intensively studied recently [Il-[1, [iTl - fl^ . It is well 
known that, for epidemic processes in scale-free networks, 
the high heterogeneity of connections leads to the ab¬ 
sence of an outbreak threshold [ itL [l8l| . In this paper, six 
different types of contact and propagation dynamics on 
a network are considered. Using degree-based mean-field 
theory, I show that the presence or absence of a threshold 
depends on the contacting and propagating dynamics. 

In this paper, I adopt the susceptible-infected- 
susceptible (SIS) model. The SIS model is one of the 


simplest models in epidemiology and is also known as 
the contact process model. In an SIS model, a popula¬ 
tion with N individuals is categorize into two compart¬ 
ments: susceptible (S) and infected (I). The disease is 
transmitted only when a susceptible individual is in con¬ 
tact with an infected individual. I assume that A is the 
rate of the contact that is enough for the disease to be 
transmitted. By rescaling time, the recovery rate can be 
set to I without the loss of generality. In the case of a 
fully mixed population, the model is represented by two 
stochastic events: 

1. At the rate of A, the contagious event is performed. 
An individual is chosen at random. If the individ¬ 
ual is susceptible, another individual is randomly 
chosen. If the second individual is infected, the 
hrst individual becomes infected. 


2. At the rate of I, the recovery event is performed. 
Choose an individual at random. If the individual 
is infected, it recovers and becomes susceptible. 


Mean field theory shows that the average density of in¬ 
fected individuals p{t) follows the rate equation 


^ = -pit) + Xpit)[l-pit)]. 

The equilibrium solution is calculated as 


P 




1-i (A>1) 
0 (A<1). 


( 1 ) 

( 2 ) 


There is an outbreak threshold Ac = 1, above which epi¬ 
demic spreads can occur. Note that if the roles of sus¬ 
ceptible and infected individuals are exchanged in the 
contagious event, result ([2]) does not change. In other 
words, the direction of the transmission is irrelevant to 
the spreading phenomena in a fully mixed population. 
In the next section, I extend the model to include the 
network structure and the direction of the transmission. 


II. MODELS AND ANALYSES 
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Here, six SIS models on networks are presented. I as¬ 
sume that an individual is located at each node of a fixed 
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network. The links of the network represent potential 
connections. It is assumed that an individual is acti¬ 
vated at each timestep. Then, the active individual can 
contact its nearest neighbors on the network. Two types 
of contacting mechanism are considered: ( 1 ) all neigh¬ 
bors or ( 2 ) only one neighbor is contacted at the same 
time. In addition, two possibilities for transmission are 
considered: an active individual is (a) the sender or (b) 
the receiver. I also consider the hybrid case (c): an active 
individual plays both roles. By combining the contact¬ 
ing and transmitting types, six models are constructed, 
as shown as follows (see also table |I| . 


To calculate the outbreak threshold, Ac, I take the limit 
0* —>■ 0, and thus, I easily obtain 


Since the second moment of the degree (fc^) diverges for 
scale-free networks with 7 < 3, the outbreak threshold 
can vanish, i.e., Ac = 0. Thus, disease can spread on 
scale-free networks no matter how low the contact rates 
may be. This model is essentially the same as the Pastor- 
Satorras and Vespignani model [ itI. [l 8 j|. 


A. Model la 


B. Model lb 


In model la, if an infected individual is activated, all 
of its neighbors become infected. The contagious event 
in the previous section is replaced as follows 

1 . At the rate of A/(fc), the contagious event is per¬ 
formed. Choose an individual at random. If it is 
infected, all of its neighbors get infected. 

Here, the contact rate A is divided by (fc) because the av¬ 
erage number of contact per unit time is set to 1. Follow¬ 
ing the degree-based mean-field theory [ itL [T 8 I|. consider 
the relative density Pk{t) of infected individuals with de¬ 
gree k. The rate equation for pk{t) is given by 

^ (3) 

The first term of the right-hand side is the recovery event. 
The second term, which represents the contagious event, 
is proportional to the contact rate X/{k), the density of 
susceptible individuals 1 — Pk{t), the degree k and the 
variable 0 (t), which is the probability that a link trans¬ 
mits the disease. Since the probability that a link points 
to a node of degree k equals kP{k)/{k), I obtain 

0(t) = ^ kP{k)pk{t). (4) 

From Eq. the equilibrium condition leads to 

,_ fcAe* 

(fc) -b kXQ* ■ ^ ’ 

The total density of the infected in the equilibrium state 
is determined as 


p*=Y,P{k)pl ( 6 ) 

k 


Substituting Eq. ([5]) into Eq. (|4]) and dividing both sides 
by 0*, I obtain 




kX 


{k) + kxe*' 


( 7 ) 


In model Ib, I reverse the direction of the contagious 
process of Model la. An active susceptible individual 
gets infected if there is at least one infected neighbor. 
The contagious event is replaced as follows 

1 . At the rate of X/{k), the contagious event is per¬ 
formed. Choose an individual at random. If it is 
susceptible and has at least one infected neighbor, 
it becomes infected. 

In this case, the transmission probability in one time step 
is 


l-(l-0(^))^ 


rather than kQ{t). Therefore, the rate equation for pk{t) 
is rewritten as 

^ + (^ [1 - [1 - (1 - ( 9 ) 

From Eq. ®, the equilibrium condition leads to 


A [1 - (1 - 0.)'=] 
{fc)+A[l-(l-0*)'=]' 


( 10 ) 


Substituting Eq. (HOI) into Eq. ®,I obtain a self consis¬ 
tent equation 


0 * 


1 

W) 


J2kP{k) 


A [1 - (1 - 0,)'=] 

(fc) + A[l-(l-0*)'=]' 


( 11 ) 


Considering the limit 0* —?> 0, it is obvious that the 
outbreak threshold is given by Eq. (jS]) once more. The 
absence of an outbreak threshold is seen as in model la. 
Since 1 — (1 — 0(t))^ < kQ{t), the equilibrium density p* 
is smaller than that for model la. 


C. Model Ic 

Next, I consider a hybrid of models la and lb. Thus, 
the contagious event is described as follows: 
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TABLE I. Properties of six models. 


model 

contacts 

active individual outbreak threshold 

equilibrium density of infected 

la 

all neighbors 

sender 

vanish 

same as Pastor-Satorras and Vespignani’s model 

lb 

all neighbors 

receiver 

vanish 

lower than model la 

Ic 

all neighbors 

hybrid 

vanish 

intermediate of la and lb 

2a 

one neighbor 

sender 

finite 

lower than model 2b 

2b 

one neighbor 

receiver 

finite 

same as well-mixed case 

2c 

one neighbor 

hybrid 

vanish 



1. At the rate of X/2{k), the contagious event is per¬ 
formed. Choose an individual at random. If it is 
infected, all of its neighbors becomes infected. If it 
is susceptible and has at least one infected neigh¬ 
bor, it becomes infected. 


In this case, the rate equation is formed as a compound 
of Eqs. dSl) and dH]): 

^ = -pUt)+^^[i-p,{t)] [kQ{t) +1 - (1 - 

( 12 ) 

In the same way as for model Ib, I obtain a self consistent 
equation 


0 * 




A [fee,-h I - (1 - 

2(/c)-f A[A:0,-f I-(l-0*)'^']' 


(13) 

Thus, the outbreak threshold is given as Eq. ([5|). The 
equilibrium density, p*, is intermediate between those for 
models la and lb. 


infected: 


pit) = ^Pik)Pk{t). 

k 

This comes from the fact that the number of transmis¬ 
sions provided by an infected individual is proportional 
to the degree in model la, while the degree is irrelevant to 
infectivity in model 2a. The equilibrium condition leads 
to 


kXp* 

{k) + kXp*' 


(15) 


Substituting Eq. (fT51) into Eq. dH]) and dividing both sides 
by p*, I obtain 


Taking the limit as p* —0, I find that the outbreak 
threshold is 


D. Model 2a 

In the above three models, it was assumed that an 
individual contacts all of its neighbors simultaneously. 
Hereafter, I introduce three models to restrict the con¬ 
tacts to only one. In model 2a, if an infected individual 
is activated, only one neighbor can become infected. The 
contagious event is replaced as follows. 

1 . At the rate of A, the contagious event is performed. 
Choose an individual at random. If the individ¬ 
ual is infected, choose another individual among 
its neighbors randomly. If the second individual is 
susceptible, it becomes infected. 

In this case, the rate equation is rewritten as 

^ = -p,(t) + A[1 - Pk{t)]^^p{t). (14) 

Here, instead of 0(t) in Eq. (I3|) for model la, I need to 
use p(t), which is the probability that an individual is 


Ac = l, 

regardless of the degree distribution P{k). 


E. Model 2b 

Here, I reverse the direction of the contagious process 
of model 2a. An active susceptible individual contacts 
only one neighbor; if the neighbor is infected, the suscep¬ 
tible individual becomes infected. The contagious event 
is replaced as follows. 

1 . At the rate of A, the contagious event is performed. 
Choose an individual at random. If the individual 
is susceptible, choose another individual among its 
neighbors randomly. If the second individual is in¬ 
fected, the first individual becomes infected. 

In this case, the rate equation for pk(t) is written as 

dpkjt) 
dt 


Pk{t) + X[l- pk{t)]Q{t). (17) 
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Here, k/{k) is removed in Eq. ((S]) for model la. The 
equilibrium condition leads to 

AP)* 

A = 

In this case, is independent of the degree, k. Thus, 
Pk — P* = ©*■ Fromeq. (ITO)) . the nonzero solution p* > 0 
satisfies 


1 + Ap* ■ 

Thus, the analytical solution is obtained as 



(19) 


The density of infected individuals p* coincides with SIS 
in a fully mixed population, regardless of the degree dis¬ 
tribution. As a result, the outbreak threshold is given 

by 


Substituting eq. (1^ into eqs. O and dH), I obtain 


0 * 

P* 




A(0* + fcpV(fc)) 


= 


2 + x{e* + kp*/{k))’ 
x{e*+ kp*/{k)) 

2 + x{e* + kp*/{k))’ 


To estimate the outbreak threshold, I take the limit as 
0 * ^ 0 and p* 0 : 


20* = A,(0*+p*(fc2)/(fc)), 

2p* = Xc{e*+p*). 

Solving Eq. (|2^ . I obtain 


l + yWJm 

For scale-free networks with exponent 7 < 3, the out¬ 
break threshold vanishes; i.e., Ac = 0. 


Ac = 1. 


III. NUMERICAL SIMULATIONS 


Jensen’s inequality leads to 

E . , kX X 

^ ^ \k)+ kXp* ~ 1 + Xp*' 

By comparing Eqs. (HH) and dl), I can deduce that in 
model 2 b, the equilibrium density p, is larger than that 
in model 2 a. 


To confirm the prediction made using the degree-based 
mean-field theory, numerical simulations are performed. 
Here, I adopt Barabasi and Albert’s model (the BA 
model) [l^ . This model assumes that a new node is added 
to the network and that each new node is connected to 
m > 1 existing nodes with a probability that is propor¬ 
tional to the degree of the existing nodes. An analytical 
approach based on the master equation shows that the 
degree distribution of the BA model is given by 


F. Model 2c 


P{k) 


2 m(m -|- 1 ) 

k{k + l)(fc -|- 2) 


(23) 


Finally, I consider a hybrid of models 2a and 2b. Thus, 
the contagious event is given as follows 

1. At the rate of A/2, the contagious event is per¬ 
formed. Choose an individual at random and then 
choose another individual among its neighbors ran¬ 
domly. If the first individual is susceptible and the 
second one is infected, the first one becomes in¬ 
fected. If the hrst individual is infected and the 
second one is susceptible, the second one becomes 
infected. 


for k > m. Thus, for large values of fc, 

P(k) ~ k-^. 

The average degree (k) is 2m asymptotically. In fig. I, 
the total density of the infected p* is plotted as a func¬ 
tion of A. The crosses represent the numerical results 
for N = 100000 and m = 2, and the curves show the 
theoretical solutions, which are calculated using Eq. (1^ 
as the degree distribution. The theoretical calculations 
agree well with the numerical results. 


In this case, the rate equation is formed as a compound 

of eqs. HZD and IHl): IV. DISCUSSION AND CONCLUSION 


dpkjt) 

dt 


-Pk{t) + ^[1- Pk{t)] 


e(i) + 


• ( 20 ) 


The equilibrium condition leads to 


X{Q*+ kp*/{k)) 

2 + x{e* + kp*/{k))' 


( 21 ) 


I have analyzed the spreading phenomena on scale-free 
networks using six SIS models with different contact and 
propagation mechanisms. The theoretical predictions 
showed good agreement with the results of the numer¬ 
ical simulations. The six models were divided into two 
types: in the first type, an active individual can contact 
all of its neighbors at the same time. Here, the outbreak 
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FIG. 1. The density of infected individnals p* is plotted as a function of A for the six different models. The curves show the 
theoretical predictions, while the crosses represent the numerical results. The theoretical results are calculated from Eqs. 0, 
([nil. (|T^. ([Tell , mi) and (1221) in the case of Eq. (I23II . In the numerical simulations, the system size is set to = 100000 and 
the average degree is set to (k) = 4. Each point is obtained by averaging over 10000 unit time after 10000 relaxation time on 
10 different network realizations. 


threshold can vanish, regardless of the direction of prop¬ 
agation. In the case where active individuals are senders 
(model la), the epidemic prevalence is larger than in the 
case where they are receivers (model lb). This is due to 
the fact that in model lb, propagation from more than 
one infected neighbors comes to nothing. In the second 
type of model, an active individual can contact only one 
neighbor. In contrast to the first type, when the active 
individuals are receivers (model 2b), the epidemic preva¬ 
lence is larger than that for model 2a. This result may 
look counterintuitive at first sight; however, this result 
is not so surprising, because for model 2b, an infected 
individual at a hub can transmit disease to more than 
one individual during a timestep, in contrast to model 
2a. It may be more surprising that in the case where 
the transmission is bidirectional (model 2c), the outbreak 
threshold can vanish, while in the case the transmission 
is one way (models 2a and 2b), the threshold remains 


finite. These results are summarized in table [H 

In conclusion, the epidemic prevalence depends en¬ 
tirely on the mechanisms of contact and propagation. In 
fact, the six models handle extreme situations and may 
not be realistic for practical contagious diseases. How¬ 
ever, the results of this paper have a wide-range of ap¬ 
plications in the study of the spreading phenomena, not 
only for epidemic diseases but also for such things as ru¬ 
mors and information. This is expected to be a reference 
point when considering more complex models, such as 
ones in which the network can change dynamically. 
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